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A long range Weak Equivalence Principle (WEP) violating force between Dark 
Matter (DM) particles, mediated by an ultralight scalar, is tightly constrained by 
galactic dynamics and large scale structure formation. We examine the implications 
of such a "dark force" for several terrestrial experiments, including Eotvos tests of 
the WEP, direct-detection DM searches, and collider studies. The presence of a dark 
force implies a non- vanishing effect in Eotvos tests that could be probed by current 
and future experiments depending on the DM model. For scalar singlet DM scenar- 
ios, a dark force of astrophysically relevant magnitude is ruled out in large regions 
of parameter space by the DM relic density and WEP constraints. WEP tests also 
imply constraints on the Higgs-exchange contributions to the spin- in dependent (SI) 
DM-nucleus direct detection cross-section. For WIMP scenarios, these considerations 
constrain Higgs-exchange contributions to the SI cross-section to be subleading com- 
pared to gauge-boson mediated contributions. In multicomponent DM scenarios, a 
dark force would preclude large shifts in the rate for Higgs decay to two photons 
associated with DM-multiplet loops that might otherwise lead to measurable devia- 
tions at the LHC or a future linear collider. The combination of observations from 
galactic dynamics, large scale structure formation, Eotvos experiments, DM-direct- 
detection experiments, and colliders can further constrain the size of new long range 
forces in the dark sector. 

o 

o> 
o 

> 

•i-H 

>< 



* seancarroll@gmail.com 

t mantry @ wise . edu 

•f mj rm@physics .wise .edu 



2 



Contents 

I. Introduction 2 

II. Fifth-Force Phenomenology 6 

III. Light Scalar Coupling to Macroscopic Objects 9 

IV. Light Scalar Coupling to the Standard Model 13 

A. Coupling to the Higgs Sector 13 

B. Non-renormalizable interactions 15 

C. Fine tuning and the light scalar mass 17 

V. WIMP DM and Eotvos Experiments 18 

VI. WEP Tests, Direct Detection, and Higgs Boson Decays 24 

A. WEP Tests and Ultralight-Scalar-Higgs Mixing 25 

B. WEP Tests and Scalar Singlet DM Relic Density 30 

C. WEP Tests and DM- Nucleus Cross-Sections 31 

D. WEP Tests and Higgs Decays 36 

VII. Dark Force Parameter Space 38 

VIII. Conclusion 39 

Acknowledgments 41 

A. Atomic Charge to Mass Ratio Under a Coupling to Ultralight Scalars 41 

B. Effective Potential for Higgs-Ultralight-Scalar Mixing 44 
References 45 



I. INTRODUCTION 

There is now compelling evidence for the ACDM model or the 'Standard Model' of cos- 
mology according to which the energy of the universe is about 74% dark energy, 22% Dark 
Matter (DM), and 4% baryonic matter. There have been independent confirmations of the 
dark energy component of the universe from observations of high redshift Type la supernovae 
[1, 2, 3, 4, 5]. The evidence for DM is even more compelling from the study of galactic rota- 
tion curves [6, 7, 8], acoustic oscillations in the cosmic microwave background [9, 10, 11, 12], 
large scale structure formation [13, 14], and gravitational lensing [15, 16]. In spite of such 
strong evidence for the existence of dark energy and DM, almost nothing is known about 
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their properties. The simplest explanation of dark energy is a small but non-zero cosmo- 
logical constant. The DM properties such as its mass, quantum numbers, and interactions 
with the Standard Model (SM) remain unknown. Furthermore, it remains to be seen if 
there is only one type of DM particle responsible for all of the observational evidence, or 
if there exists a rich spectrum of DM particles analogous to the complexity seen in the 
visible sector. Many experiments are underway to detect DM and determine its proper- 
ties. Ground based direct detection experiments [17, 18] put limits on the DM mass and 
the strength of its interaction with baryonic matter from observations of recoiling nuclei. 
Experiments [19, 20, 21, 22] studying cosmic rays from the galactic halo have recently seen 
indications of an electron/positron excess, which could be interpreted as evidence for DM 
annihilation, and can constrain the DM mass and interactions. There has also been a re- 
cent proposal to observe a possible DM magnetic moment via the gyromagnetic Faraday 
effect [23, 24]. 

Another set of experiments are devoted to question of whether DM violates the Weak 
Equivalence Principle (WEP) and are the focus of this paper. There exist a variety of 
scenarios for new interactions confined solely to the dark sector and the possibility that 
they might lead to WEP violation. The possibility of gauge or Yukawa forces confined to 
the dark sector have been studied in other contexts [25, 26, 27, 28, 29, 30, 31, 32]. Short 
range WEP tests have been studied in [33, 34, 35] for example. In this work, we focus on 
a long range dark force, mediated by an ultralight scalar, and study its implications for 
terrestrial experiments. For this scenario, the dark force can be communicated to ordinary 
matter via virtual DM loops that connect the scalar with ordinary matter, as long as the DM 
candidate is not sterile. This mechanism will give rise to effects in terrestrial experiments. 
We investigate the resulting impact that WEP violation in the dark sector may have on 
DM detection experiments, laboratory based WEP tests, or even studies of Higgs boson 
properties at colliders. Constraints on WEP violation in ordinary matter induced by WEP 
violation in the dark sector were recently studied in [36, 37]. In addition, a connection 
between direct DM detection experiments and WEP violation was shown in [36]. 

Many models that contain the interaction of an ultralight scalar with DM [38, 39, 40, 
41, 42, 43, 44, 45, 46, 47, 48, 49] have been proposed to explain features in the DM dis- 
tribution and explore the possibility of DM-quintessence interactions. More recently, work 
with non-universal scalar-tensor theories of gravity with the Abnormally Weighting Energy 
(AWE) Hypothesis [50, 51] also invoke couplings of an ultratight scalar to the dark sector as 
a way of explaining the observed cosmic acceleration even in the absence of a dark energy 
fluid. Constraints on such scenarios from big bang nucleosynthesis have also been stud- 
ied [52] . There are several other observational motivations, including higher than predicted 
supercluster densities [53] and voids [44, 48] ( for a summary see [36, 54]). The existence 
of a long-range attractive Yukawa force between DM particles would accelerate structure 
formation and could help explain some of these observations. Strong constraints on such a 
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dark force are derived from observations of DM dynamics in the tidal stream of the Sagit- 
tarius dwarf galaxy [54, 55], which indicate a force with strength less that 20% of gravity. 
However, new observational systematic errors have been recently discovered [56] that could 
require a revision of this result, perhaps allowing for a stronger dark force. A more recent 
analysis [57] considers the effect of a dark force on the evolution of density perturbations 
and the resulting impact on the CMB spectrum. This analysis constrains the strength of a 
dark force to be less than 5% of gravity. 

From a purely theoretical perspective, the existence of an ultralight scalar <fi with mass 
m<p < 1CT 25 eV, able to mediate a long-range force over scales of interest to galactic dynamics, 
would introduce a new hierarchy in addition to that between the weak scale mw ~ 100 GeV 
and the Planck scale M P ~ 10 19 GeV. However, as we still await experimental evidence 
for a mechanism to explain the the hierarchy between the weak and Planck scales, and in 
light of the discovery of an unnaturally small cosmo logical constant, we keep an open mind 
and do not attempt to provide an explanation for the ultralight scalar mass. We assume 
the existence of a finely-tuned ultralight scalar mediating a long range force dark force and 
study its consequences for terrestrial experiments. 

In what follows, we amplify on our earlier work[37] and that of Ref. [36], using simple 
DM scenarios to illustrate the prospective implications of dark sector WEP violation for 
terrestrial experiments. We study three representative minimal DM scenarios to explore 
the range of possible implications: scalar DM that is a singlet with respect to SM gauge 
interactions; scalar DM that is the neutral component of a real SU(2)x, triplet; and fermionic 
DM that lives in a vector-like representation of SU(2) L . Our main conclusions are: 

(i) The presence of a dark force implies a non-zero effect in Eotvos experiments if the 
DM interacts with Standard Model (SM) fields. For scalar singlet DM, this effect 
arises from loop-induced mixing between the ultralight scalar and the SM Higgs, while 
for representative WIMP scenarios (scalar or fermionic) additional contributions arise 
from loop-induced higher-dimension operators that couple the ultralight scalar directly 
to matter. We derive or der-of- magnitude expectations for the minimum size of this 
effect for these representative scenarios as illustrated in Fig. 4 for WIMP DM and 
Table. II for scalar singlet DM. For a dark force with strength roughly 20% of gravity, 
one could expect a non-vanishing effect, for non-minimal WIMP DM models and 
in certain regions of parameter space of scalar singlet DM models, within reach of 
future approved Eotvos experiments such as Microscope [58] able to detect anomalous 
accelerations to a sensitivity of Aa/a ~ 10~ 15 . The MiniSTEP experiment [59] with 
an increased sensitivity of Aa/a ~ 10~ 18 , currently under study by NASA and the 
ESA, could see non-vanishing effects in minimal WIMP models which can induce a 
WEP violating starting at two loops. 

(ii) For scalar singlet DM, a dark force of astrophysical relevance, is already ruled out in 
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large regions of parameter space. The bounds from Eotvos experiments constrain the 
size of DM-Higgs interaction which determines the relic density. In large regions of 
parameter space, the bound on DM-Higgs interactions implies a suppression in the 
DM annihilation rate resulting in a relic density that over-closes the universe. As a 
result, relic density considerations in scalar singlet DM models can yield the strongest 
bounds on the size of a dark force. 

(iii) The constraints on DM-Higgs interactions lead to upper bounds on the magnitude 
of Higgs-exchange contributions to DM-nucleus cross sections. These bounds depend 
on the Higgs mass, implying that a combination of direct detection experiments and 
Higgs boson discovery could be used to test the simplest scenario for a dark force 
in some DM scenarios. In particular, Higgs-exchange contributions dominate the SI 
scalar singlet DM-nucleus cross section, so that dark force considerations - together 
with the observed DM relic density - imply constraints on the entire cross section. In 
contrast, WIMP-nucleus cross sections receive contributions from electroweak gauge 
boson-exchange that are not constrained by the presence of a dark force. As we show 
below, dark force considerations and present limits from Eotvos experiments imply 
that the Higgs exchange contributions are sub-leading compared to those from gauge 
boson-exchange. The corresponding bounds for the scalar singlet and real triplet 
models are illustrated in Figs. 10. Tests of the WEP can only constrain the full DM- 
nucleus cross section if the DM particles are singlets with respect to the SM gauge 
symmetries (see, e.g., [60, 61, 62] and references therein) so that elastic scattering 
proceeds only via Higgs exchange (at least at tree level). 

(iv) In multi-component WIMP DM scenarios, where one of the light (^ 200 GeV) DM 
components has a non-zero coupling to the Higgs, the presence of a dark force - 
together with tests of the WEP - imply testable upper bounds on one-loop WIMP- 
induced shifts in the branching ratio for the SM Higgs to decay to two photons. 
These bounds generally lie well below the prospective sensitivities of LHC studies 
of Br(H — > 77) as seen in Fig. 12. The observation of a significant shift in this 
branching ratio would likely preclude this scenario for a dark force. 

(v) The existence of an observable long range dark force which requires < 10~ 25 eV, 
implies restrictions in the space of finite renormalized parameters in addition to the 
usual fine tuning of radiative corrections that are sensitive to the cutoff. We discuss 
these regions in parameter space and their implications for the observation of a dark 
force. 

In arriving at these conclusions, we emphasize we have drawn upon representative cases 
rather than carrying out a comprehensive study. We expect that our conclusions will gen- 
eralize to other DM scenarios, but do not preclude the possibility of exceptions in some 
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cases. We also note that our analysis and conclusions differ from those of Ref. [36], who 
first observed that bounds on WEP and the presence of an astrophysically relevant dark 
force could imply constraints on DM-nucleus cross sections. The bounds obtained in that 
work lie well below the reach of future direct detection experiments. In what follows, we 
argue that an effective operator analysis consistent with the fine-tuning needed to maintain 
a vanishingly small scalar mass implies considerably weaker bounds than given in Ref. [36] . 

The outline of the paper is as follows. In section II we review the phenomenology of 
experimental WEP tests and establish notation. In section III we review the derivation of 
the ultralight scalar coupling to macroscopic objects in terms of its couplings to the Standard 
Model (SM) particles. In section IV we discuss in a model independent way the mechanisms 
by which the ultralight scalar can couple to the SM. In sections V and VI we examine the 
experimental consequences of a dark force for various minimal DM models. In section VII 
we discuss the regions in parameter space where an observable dark force is possible and 
how they relate to our analysis. We conclude in section VIII. 



II. FIFTH-FORCE PHENOMENOLOGY 



We begin by considering the force between two bodies mediated by a scalar field with 
mass m^. In the non-relativistic limit, the Yukawa potential between a test body i and a 
source s separated by a distance r is given (in units where % = c = 1) by 

V« = -tes^e-^, (1) 

where Q i>s denote the charges of the test and source objects under the force mediated by <fi. 
The parameters £ ijS are 1 

{1 for fermionic objects, 
2^— for scalar objects. 

Note that the charges Qi yS are of mass dimension one and zero for scalar and fermionic objects 
respectively, so that the equation is dimensionally consistent. These mass dimensions will 
become apparent when we study specific models. The Newtonian gravitational potential 
between a body with mass Mj and a source with mass M s is 

V„ = -™™, (3) 



1 The t-channcl <p exchange amplitude is accompanied by an extra factor of 2rrii, s for fermions relative to 
scalars. This is due to the fermionic spinor normalization Ui yS Ui yS = 1m i s in the non-relativistic limit. 
These factors are absorbed by switching to states with normalization (p|q) = (27r) 3 5 < ' 3 - ) (p — q) in order to 
compare with the non-relativistic Born amplitude. For scalars we are then left with an additional factor 
of in the potential relative to fermions. 
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where G is Newton's constant. It is therefore convenient to write the total potential as 

V = - GM ^ {1 + a is e-™* r ) , (4) 

where 

oti a = ^—^ Kits, (5) 

47TG [lifis 

is a dimensionless parameter characterizing the strength of the new force relative to gravity, 
expressed in terms of the charge-to-mass ratio qj\i = Q/M, where \i is the mass in atomic 
mass units. The parameters £ is are 



for fermionic objects, 
for scalar objects. 



(6) 



The parameter oti S is not universal and in general depend on the composition of the macro- 
scopic bodies acting as sources for <fi. 

Eotvos experiments look for violations of the equivalence principle by measuring the 
difference in acceleration of two test bodies of different compositions in the presence of a 
common source. Experimental constraints on new long-range composition-dependent forces 
are typically expressed in terms of the Eotvos parameter, 

Aa 



V = 2 



Oi — d2\ 



(7) 



\ai + a 2 \ 

where Oj is the total acceleration of object i — 1,2, Aa = a\ — a 2 , and a is the universal grav- 
itational acceleration in the absence of any new long range forces. The last approximation 
made above is valid when the fifth force is weaker than gravity. From (4), the acceleration 
of object i due to the source s is 



GM„ 



[I + a is {l + m r)e-^ r ] . 



(8) 



We are interested in forces that are considerably weaker than gravity, and distances less 
than the Compton wavelength of the scalar, r mT 1 . The Eotvos parameter is then 



vl' 2 



AuG 



(9) 



Currently, the strongest limits on violations of the weak equivalence come from torsion 
balance Eotvos experiments [63] which give the constraints 



Be,Ti 
'E 



< (0.3 ±1.8) x 10 



-13 



Be,Ti 

'DM 



< (4± 7) x 10~ 5 . 



(10) 



The Eotvos parameters r/^ e > Tl and ?7^' Tl measure differential acceleration of laboratory test 
samples of Beryllium and Titanium with the Earth and galactic dark matter as the source 
bodies respectively. 



8 



Experiment 


Expected Future Sensitivity in i] 


MiniSTEP[59] 


10" 


-18 


Microscope [58] 


10" 


-15 


Apollo (LLR)[64] 


10" 


-14 



TABLE I: Expected sensitivities for the Eotvos parameter in future experiments testing the WEP. 
The MiniSTEP experiment is currently under study by NASA and the ESA. Microscope has been 
approved and the Apollo(LLR) experiment is underway. 



Future experiments, currently being studied, are expected to further improve the the 
bound on the Eotvos parameter by several orders of magnitude as shown in Table. I The 
MiniSTEP experiment [59], currently under study, would use test objects of different com- 
position orbiting earth in free fall and new technology to reduce thermal noise. If approved, 
this experiment is expected to achieve the highest sensitivity of rj ~ 10~ 18 . The Microscope 
experiment, which has been approved, uses the same principle but is expected to reach a sen- 
sitivity of rj ~ 10~ 15 . In the method of Lunar Laser Ranging (LLR) used by the APOLLO 
collaboration [64], the differential acceleration of the Earth and Moon is measured in the 
presence of a source like the Sun or galatic dark matter. The APOLLO collaboration, which 
is currently underway, is expected to achieve a sensitivity of rj ~ 10~ 14 improving the current 
bound on rj by an order of magnitude. Methods using atom interferometry [65] could reach 
a sensitivity of rj ~ 1CT 17 . 

These experiments are also sensitive to WEP violation in the dark sector, if the DM has 
interactions with the SM. Through quantum effects involving virtual DM, WEP violation 
in the dark sector will be communicated to ordinary matter and these effects can be tested 
in Eotvos experiments. WEP violation in the dark sector is already constrained from an 
analysis of the tidal disruption in satellite galaxies [54]. This study constrains the coupling 
of to DM particles by putting bounds on the parameter (3 



where g x denotes the DM charge under the fifth force, M x denotes the DM mass, M P = 
1/VG is the Planck mass, and £ x is as denned in Eq. (2). The coupling g x appears in the 
Lagrangian via interaction terms for fermionic 2 and scalar DM of the form 



= J 9 x XX<t>, fermionic DM, 
9xX ] X<t>, scalar DM, 



Thus, we see that for fermionic DM, g x is dimensionless and for scalar DM it has dimen- 



2 For simplicity we assume that the fermionic DM is in a vector-like gauge representation so that xx4> 
is gauge invariant. In the more general case the coupling of <f> to fermionic DM may arise from higher 
dimension operators. 



9 



sion one. From the analysis of tidal streams in the Saggiatrius galaxy, Kamionkowski and 
Kesden[54] obtained the approximate upper bound of 

13 < 0.2. (13) 

Newly discovered systematic errors [56] could lead to a revision of this bound and more 
recently, the work of [66] showed the possibility of f3 ~ 1 consistent with observations 
of galactic dynamics. A more recent analysis [57] of the CMB and large scale structure 
formation gives a tighter bound of f3 < 0.05. In this paper we use (3 = 0.2 as a reference 
value for most discussions, and our results be straightforwardly translated to other values 
of (3. 

III. LIGHT SCALAR COUPLING TO MACROSCOPIC OBJECTS 

The charge to mass ratio under a fifth force for an elementary particle is straightforward 
to obtain in terms of the Lagrangian parameters. For example, the charge to mass ratio for 
elementary fermionic or scalar DM x is given by 

(-) (14) 

This charge to mass ratio is obtained by computing the tree level exchange diagram 
between two DM particles and taking the non-relativistic limit to compare with Eq. (4). 

For composite materials the calculation of the charge to mass ratio is more complicated, 
as one has to take into account hadronic, nuclear, and atomic matrix elements of various 
operators containing SM fields that couple to <fi as well as the effects of binding energy. In 
particular, one needs the charge to mass ratio for the various types of atoms that make up 
the laboratory test materials. We compute these ratios using an effective field theory valid 
near the nucleon mass scale that involves the light quarks q = {u, d, s}, gluons, the charged 
leptons £ = {e,/i}, the photon, and the light scalar <p. All other heavier degrees of freedom 
have been integrated out. The interaction terms in this effective Lagrangian take the form: 




As we discuss below, the effects of the coupling to heavy quarks, the tau lepton, massive 
gauge bosons, and x that have been integrated out are encoded in the operator coefficients 
g q / and c gn . We assume that the couplings of to the SM fermions are linearly proportional 
to the fermion mass. This will make the analysis simpler, as we will see, by allowing us to 
exploit the scale invariance of the energy momentum tensor. This assumption is realized in 
several types of DM models. The couplings c g and c 7 can be straightforwardly computed in 
any given model. To illustrate, consider a model in which <p couples to the SM fermions at 
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the electroweak scale as 

£<t>ff = — m q QQ4> + — m e Zty + —m T TT<f> + —rriQ QQ4>, (16) 
m m m m 

q I Q 

where the sum over Q denotes a sum over the heavy b, c, t quarks and all the couplings 
9e,r,q,Q above are independent of the SM fermion masses. One can then integrate out the 
heavy quarks to obtain [67] the renormlaization group invariant relation 

fy fy ry 

m Q QQ = ~G a ^ ~ 1^ V". rnrfr = V" (17) 

to leading order in the heavy quark expansion and perturbation theory. Note that the RHS 
above is independent of the heavy quark mass. In this case the couplings c g and c 7 in 
Eq. (15) are given by 

at leading order. The mass operators on the LHS of the equations in Eq. (17) appear in the 
QCD+QED energy momentum tensor and are scale invariant, allowing us to evaluate a s 
and a in Eq. (18) at the low energy scale of the effective theory (when taking the nucleon 
matrix element). Because they do not run below the electroweak scale, couplings ^q ;T are 
evaluated at that scale. 

We now evaluate the coupling of <fi to an atom [35] of type 'A'. Doing this allows us to 
determine the charge to mass ratio qA^Al Ha needed for Eotvos parameters, as seen from 
Eq.(9), if the test or source bodies are made up of atoms of type 'A'. We define the effective 
atomic coupling as 

^ f g A AA(f), fermionic atoms, 

I gAA^A(fi, scalar atoms, 

where the A is the field that destroys the atomic state and again qa is dimensionless for a 
spin 1/2 atom and has dimension one for a spin zero atom. We determine qa by a matching 
calculation 

(A\£ A aM<P) = 9 A U = {AlCMti, (20) 

where we have used a non-relativistic normalization for the atomic states (A(p)\A(q)) = 
(2n) 3 5 3 (p — q), £a is the normalization factor defined in Eq. (2), and is defined in (15). 
From Eq.(20), as explained in Appendix A, the general expression for the charge to mass 
ratio q A iAl^A is 




9aU = 2c g g 3 J_ 
M A ' p 3 M A 



Z(( e m e + ^ Q m q x q,p) + ( A ~ Z ) XI x ^ n + Ua 



(21) 
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where the quantity uja is given by 



u A = k{A\F^F, v \A) - Y.^ mk 



dS A 
drrik ' 



(22) 



£a is the atomic binding energy as defined in Eq.(A5), the quantities ( k and k are given by 



G 



9k 2# 3 



m r , 



K 



as in Eq.(AlO), and x q>p and x q;Tl denote the nucleon matrix elements 

x q , P = {p\qq\p}, x q ,n = (n\qq\n) , 



(23) 



(24) 



which are known experimentally [68, 69] and given in Eq.(A13) of appendix A. In Eq.(23), 
(3 3 and f3 e denote the QCD and QED beta functions respectively. 

Using Eq.(21) in Eq.(9) for test objects made up of atoms with atomic weights A x and 
A 2 , the general expression for the Eotvos parameter r] s with source S is 



V s = 



Ml 



[lJ s 



■ Z 1 

'Ma, 



M A . 



-) (( e m e + (q m q x q , p ) 



+ ( 



A 1 -Z 1 A 2 - Z 2 



M A 



Ma 2 



M Al M A2 



(25) 



where {^) s denotes the charge to mass ratio for the source object and A k , Z k (k = 1,2) 
refer to the atomic weights and atomic numbers of the two laboratory samples. For order of 
magnitude estimates, we follow Ref. [35] and ignore binding energy effects, encoded in the 
quantitities uja 1 , 2 - Setting Ma — ArriN for the atomic masses, we then obtain the simpler 
expression 



M 2 P 



N 



fi/ s 



^ ~ T 2 ^^ eme + zZ^ m i ( X «>P ~ X 9,n)} 



(26) 



From Eqs. (18) and (23), the parameters ( k appearing above are given by 



Ck 



1 



9k 



Q 



(27) 



at leading order. Here g k denotes the couplings of <p to the light quarks and charged lep- 
tons and gQ denotes its coupling to the heavy (b, c, t) quarks. A special case that will be 
of particular interest in subsequent discussion occurs when the couplings to fermions are 
universal, apart from the fermion Yukawa couplings explicitly factored out via the factors 
of rrif in Eqs. (15,16). Setting 

9k = 9q=9 (28) 
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and m p = m n = tun leads to 



~ Q ( 



\47rm 



jJL/ S 



(j~ ~ J^) { m e + Yl m i ^ ~ X ^ } 
■ ! 2 q 



(29) 



Typical source objects £ 5" used in Eotvos experiments include the Earth, the Sun, and 



galactic DM, and one needs to obtain their charge to mass ratio £c 



that appears in 



Eq.(29). If galactic DM is made of of elementary particles, then as already discussed, the 
charge to mass ratio under the dark force is given by 
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Q 



fJ,/ s 



9x_ 
M v 



(30) 



S=DM 

For objects like the Earth that are made up of many different types of atoms, the effective 
charge-to-mass ratio is obtained by a superposition of the couplings of <fi to all the different 
atoms present in the object. In contrast to the situation for differences in charge-to-mass 
ratios for test bodies, it suffices to approximate this ratio for the bulk source object by 
ignoring atomic binding energy effects and summing over the couplings of <p to all the 
neutrons, protons, and electrons present. Doing so in the case of the Earth leads to 

9pNp + 9nN n + g e (m e /m N )N e 
.fi/E m N (N p + N n ) + m e N e 

where g p and g n denote the couplings of <f> to protons and neutrons, respectively: 

g N = (N\£„j\N) , 

for iV = p or n. In the limit of a universal coupling as in Eq. (28), we have 

v 



Vu/ , 



(31) 



(32) 



9n = 9hg 



m N 



(33) 



where v = 246 GeV is the vacuum expectation value of the neutral component of the Higgs 
doublet and gh/V2 is the coupling of the physical Higgs boson to the nucleon 



V q Q / 



(34) 



Using similar methods to those employed to determine gA and ignoring small difference 
between the neutron and proton coupling, one has [70, 71, 72] 



g h ~ 1.71 x 10- 3 . 
The resulting expression for the Earth's charge-to-mass ratio in this case is 



E 



Q 



[IJ E 
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g h (N p + N n ) + {m e /v)N e 
m 2 N J (N p + N n ) + (m e /m N )N e 



0.0017^ 



m 



N 



(35) 



(36) 



The number of protons, neutrons, and electrons are N p ~ 1.9 x 10 51 , N n ~ 2.0 x 10 51 , and 
JV e ~ 1.9 x 10 51 respectively. We will make use of Eq. (36) in what follows. 
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IV. LIGHT SCALAR COUPLING TO THE STANDARD MODEL 

We now give a model-independent discussion of the coupling of ordinary matter to the 
ultralight singlet scalar that mediates the long range force. In doing so, we will lay the 
groundwork for calculating the parameters gj (/ = q,Q,£), c g , and c 7 of Eq. (15) and 
Eq.(16) or equivalently the parameters (k and k in Eqs. (23) and (27). In general there 
are two mechanisms for a singlet scalar to couple to the SM fermions and gauge bosons. 
The first mechanism involves mixing between the ultralight scalar and the Higgs, which 
allows the ultralight scalar to couple to the SM fermions and gauge bosons. The second 
mechanism entails a coupling of the ultralight scalar to the SM through higher-dimension 
(non-renormalizable) operators. We discuss these two mechanisms in this section and estab- 
lish notation. We also address the need for fine-tuning of the ultralight scalar mass when 
its interactions with the SM are non-negligible, looking ahead to a similar issue when we 
consider its coupling to DM. 

A. Coupling to the Higgs Sector 

We assume that the mediator of the dark force carries no SM charges and that it can be 
described by a gauge singlet S. There exist no renormalizable couplings of such a singlet 
scalar to the SM fermions or gauge bosons, but it can couple to the SM Higgs doublet 
with operators of mass dimension n < 4. After electroweak symmetry breaking, the n = 3 
interaction H^HS will generate mixing between S and the neutral component of the Higgs 
doublet, h. We will identify the ultralight force-carrying scalar with the lighter mass 
eigenstate, and the heavier eigenstate with the physical Higgs boson. The Lagrangian for 
the singlet S including its renormalizable and super-renormalizable interactions is given by 

C=±d ll SdvS-V(H,S), (37) 

where the potential is 

V(H,S) = -fi^H + j(H*H) 2 + S -±H*HS + ^H*HS 2 

- (^)s+|s 2 + |s 3 +^? 4 . (38) 

We have shifted the scalar S so that it has no tree level vacuum expectation value without 
loss of generality. We follow the notation of Ref. [60, 73], which explored the presence of 
such a singlet scalar in the context of collider phenomenology. The parameters may arise 
from a more fundamental theory of which the S is a residual, low-energy degree of freedom. 
As we discuss below, they may also receive contributions from DM loops if the DM particles 
couple to the both H and S. 
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FIG. 1: Interaction of S with SM fermions by mixing with the Higgs via the operator SH^H. Here, 
"X" denotes the vacuum expectation value of the Higgs doublet. 



After electroweak symmetry breaking the WHS interaction induces mixing between the 
Higgs boson h and the scalar S. In unitary gauge the neutral component of the Higgs 
doublet H is given by 



H 



o _ v + h „. _ . IWh 



V2 



and the mass terms in the potential are 

1 



Knass = ^ h h + n s S + (i hS hS), 



where 



2 2 

= — > V>S = K 2 + 



/4s = ^ 



(39) 



(40) 



(41) 



The mass eigenstates h± in terms of S and h can be written in terms of a mixing angle 9 as 

.r 



h- = Scos9 — /isin#, fo + = Ssin^ + hcos9, tan# 



1 + Vl + x 2 ' 



(42) 



with corresponding masses 



2 _ /4+j4 , /4^/4 

flL_L_ — Zl 



and we have denned 



/4s 



Z 1 ^ A*l 



(43) 



(44) 



We assume that m 2 <C m5j_, so that the physical Higgs boson and light scalar are h + and 
/i_ respectively. The light scalar h- can couple to quarks and charged leptons through its 
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mixing with the Higgs as shown in Eq. (42). We identify the light scalar that mediates the 
long range force as 

= h-, = m_. (45) 

The scalar couples to the SM fermions through its Higgs component, giving rise to the 
couplings <?/, where / denotes any of the light quarks q — u,d, s, charged leptons £ — e, fj,, r, 
or heavy quarks Q — c,b, t. One has 

g f = -sin#^^ = -sin#^, 46 

rrif v v 

where the extra factor of m p /mf after the first equality is included to be consistent with the 
convention in Eq. (15). This process is depicted in Fig. 1. We see that in this mechanism the 
coupling of (j) to ordinary matter is proportional to sin 9, with the constant of proportionality 
given entirely in terms of known quantities. The mixing angle 9 will also receive corrections 
at the loop level and in the rest of the analysis we assume that 9 is the renormalized mixing 
angle. 

For later use, we note that in the limit that fj,h — rah ^> //# corresponding to a small 
mixing angle 9, we can write 

<^s-^- (47) 

The existence of an ultralight scalar that can mediate a dark force over intergalactic distances 
requires < 1CT 25 eV. In addition to the usual fine tuning of the parameters us and fj,hs 
against radiative corrections sensitive to the cutoff (see section IV C), the finite renormalized 
parameters /is and /ihs are restricted in parameter space to satisfy the condition < 10~ 25 
eV in Eq.(47). As we will discuss in section VII in more detail, this gives rise to three types 
of regions in parameter space. In the first region, ji s and /ihs are both individually small in 
which case there will be no observable dark force. In the second region, /is and /ihs are large 
enough to give rise to an observable dark force but cancel against each other in Eq.(47) to 
maintain an ultralight mass. In the third region, as will become clear in later sections, us 
and HhS are again individually small as in the first region, but each is determined by a sum of 
much larger terms that cancel among each other. The second region is phenomenologically 
the most interesting and is the focus of this paper. 



B. Non-renormalizable interactions 



If both the S and H couple to additional fields with masses above the electroweak scale, 
then these interactions will in general induce higher dimensional operators that involve both 
the S and H in a low-energy effective theory that does not contain the heavy degrees of 
freedom explicitly. Minimal dark matter models, for example, can require TeV-scale DM 
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particles in order to achieve the observed relic density [74], and these fields may generate the 
higher dimensional S-H operators. At dimension five, one has seven independent operators 
coupling S to the SM fermions and gauge bosons: 

Of = SQ L eH^C^u R + h.c, 

Of = SQ L HC*d R +h.c, 

Of = SL L HC?e R +h.c, 

w = c w STr[W^W^], 

q b = C B S B^B^ U , 

O g = C G S Tr [G^G^], 

h = c H S{H^H)(H^H). (48) 

The flavor indices on the fields Ql, L,u R ,e R and the matrices Cf de are suppressed for 
simplicitly. Operators of the form S Ql^PQl, S u R ipu R , S d R ipd R , and S e R ipe R can be 
related to the operators 0^ d by using the equations of motion 

ipQ L = eH^Y u u R + HY d d R , 
ipu R = eHYjQ L , 
ipd R = H^YjQ L , 

ipe R = &Y}L L , (49) 

where Y) denotes the matrix of SM Yukawa couplings. We have omitted operators that 
involve derivative or pseudoscalar couplings of S. Such couplings are spin dependent and 
have a negligible effect in experiments which use unpolarized test objects. 

In general, the Wilson coefficients Cf d e are 3x3 matrices in flavor space, and can lead 
to flavor changing interactions of quarks and leptons with S. Since the couplings of S to 
quarks and leptons are extremely small (as dictated by the WEP violation bounds) there is 
no danger of introducing dangerous flavor-changing neutral currents. In the specific model 
examples considered in subsequent sections of the paper, find that the Cf d e are proportional 
to the Yukawa matrices: 

Cu,d,e = c u,d,eY u ,d,e, (50) 

where c U) d, e are the constants of proportionality. After expressing the fermion fields in the 
mass basis, in unitary gauge where the operators 0^ d become flavor diagonal. We can 
write 

6f = C uV \S u\ H° u R + h.c = cj, S u l L H° U R + h.c, 
Of = c d y d S d\ H° d R + h.c ee c\S d\ H° d R + h.c, 
Of = c e ylS 4 H° e R + h.c = & e S? L H° e R + h.c, (51) 
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where H° is the lower component of the Higgs field H in unitary gauge before electroweak 
symmetry breaking, the index % = {1,2,3} runs over the three flavor generations, and we 
have defined c l a = c a y l a . 

After electroweak symmetry breaking, the operator O h in Eq. (48) will also contribute 
to sin 9. To linear order in h one has 

O h -> C H v 3 Sh , (52) 

thereby generating a contribution to the off-diagonal element of the mass-squared matrix 

A^ hS = 2C H v 3 . (53) 

We will explore the consequences of this term when discussing scalar DM models below. For 
the moment we assume that this contribution has been included in sin 6*. 

Collecting the contributions to the couplings g q ,e,Q from the higher-dimension operators 
0^ d e and mixing effects after electroweak symmetry breaking, the coupling of the ultralight 
scalar to SM fields at the electroweak scale is given by 



m 



771/ 



9f( v ) = — ~ cos6^-=c\(v ) — sin6> — ~{y) 



V2 



~ ^^-c f (v)- sinewy), (54) 

771/ y/2 7 V ^ ' 

where the last approximation is obtained from cos 9 ~ 1 since 9 is constrained to be very 
small. We have included an extra factor of m p /mf on the RHS above to be consistent 
with the convention in Eq. (15). We have ignored the running between the scales A ~ TeV 
and the electroweak scale for simplicity, but these effects can be incorporated by computing 
the appropriate anomalous dimension matrix and solving the corresponding renormalization 
group equations. We can now use Eq. (54) in Eqs. (23) and (25) to compute the Eotvos 
parameters. In particular, we note that the contribution proportional to sin# is universal, 
so its contribution to the rjs can be evaluated using Eqs. (29) and (36). 

In general, the origin and parametric dependence of sin 6? and c/ are independent. In most 
of the parameter space where there are no strong cancellations between the two terms in 
Eq.(54), WEP violation constraints can separately bound each of the two terms in Eq. (54). 
In the next three sections we use this feature with the representative minimal DM models, 
in the presence of a dark force mediated by 0, and determine the implications for terrestrial 
experiments of direct DM-detection, Eotvos experiments, and the colliders. 



C. Fine tuning and the light scalar mass 



Before proceeding, we observe that in the absence of a symmetry that protects the light 
scalar mass from significant renormalization, one must resort to fine tuning to maintain the 
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long range character of the dark force. To illustrate, we consider the contributions from the 
WHS and H^HS 2 in V(H, S) to the singlet self energy. In the unitary gauge one has 



5 2 

£(p 2 )fftff S = 



1287T 2 



- - 7 + ln47T + ln/i 2 - F (m 2 H ,m 2 H ,p 2 ) 

e 



- 7 + 1 + In 47T - In 



(55) 
(56) 



where we work in d = 4 — 2e dimensions, \i is the corresponding t'Hooft scale, and 

F (a,b, c) = / dx In [(1 — x)a + xb — x(l — x)c] . (57) 
Jo 

Had we regulated the integrals with a momentum cutoff A co , the quadratic divergence 
proportional to 52m 2 H /e would be replaced by an expression proportional to (^A^o- 

For either choice of regulator, preservation of a tiny scalar mass requires a mass coun- 
terterm 5fi 2 s to cancel both the quadratic and logarithmically divergent contributions. In 
addition, as one sees using dimensional regularization, the finite, /^-dependent contributions 
require a corresponding /x-dependence in <5/x| as needed to maintain the scale-independence 
of the physical pole mass (or range of the dark force). The divergent and /x-dependent finite 
contributions can be minimized by either choosing 8± t 2 to be sufficiently tiny or by allowing 
for large cancellations between and the one-loop contributions. For the particular 

example discussed here, the finite contributions can also be minimized by choosing /i m m#, 
but these contributions will not be small at all scales unless the coefficients 5± t 2 are tiny or 
there exists a large cancellation (fine-tuning) between <5/z|(/i) and the one- loop contributions 
for /i ^ m H . In short, allowing for any appreciable interaction between the singlet S and 
the Higgs sector of the SM invariably requires fine-tuning at some scale in order to ensure 
that the dark force mediator remains ultralight. In what follows, we will return to this 
point when considering the coupling of S to DM. Even after allowing fine tuning, the finite 
renormalized parameters are restricted in parameter space in order to maintain an ultralight 
mass for the dark force mediator. We will discuss this issue in more detail in section VII. 



V. WIMP DM AND EOTVOS EXPERIMENTS 



In an earlier work [37], we examined constraints on the size of the coupling of an ul- 
tralight scalar to ordinary matter induced via virtual WIMP DM. The connection between 
constraints from galactic dynamics and Eotvos experiments and the size of the ultralight 
scalar couplings to DM and ordinary matter were shown in Fig. 1 of [37]. In analyzing 
the astrophysical constraints, we assumed only an upper bound /3 < 0.2 and showed that 
in representative WIMP scenarios, it leads to stronger constraints on the strength of the 
0-WIMP coupling than do the present Eotvos bounds on rj EDM . An improvement [37] of 
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about eight orders of magnitude in Eotvos experiments would be required to compete with 
the bounds from astrophysical constraints. 

Here, we explore the prospective implications of a non-vanishing (3. The presence of a 
modified, long-range dark force could help alleviate tensions in the ACDM paradigm (we 
refer the reader to Refs. [36, 54] for an extensive discussion). In what follows, we show that 
a non- vanishing (3 implies a lower bound on r) E DM in simple WIMP scenarios, so that future 
Eotvos experiments with improved sensitivity could be used to test this possibility. For 
purposes of illustration, we consider both scalar and fermionic WIMP DM. For fermionic 
WIMPs we restrict our attention to vector-like gauge representations, which simplifies the 
structure of the coupling of the ultralight singlet scalar to DM. 

The Lagrangian for minimal WIMP DM takes the form 



where c = 1/2 for a real scalar and c = 1 for a complex scalar. The covariant derivative 
depends on the SU(2) L and U(l)y representations of x- Assuming that a single WIMP 
species saturates the relic density, one finds that typical masses of such DM candidates are 
in the TeV range [74]. In general, V(x,H) can contribute to the scalar DM mass after 
electroweak symmetry breaking. However, since the typical WIMP DM masses are in the 
TeV range, such a contribution will be much smaller than the size of the mass parameter 
M ~ TeV in the second line of Eq.(58). In what follows we assume this parameter M to 
be the total DM mass since we are only interested in order of magnitude estimates. For 
gauge singlet scalar DM models with DM masses in the 100 GeV range, the contribution 
to the mass from electroweak symmetry breaking can be important. We will consider the 
case of singlet scalar DM in the next section. Furthermore, since electroweak symmetry 
breaking can in general induce mixing between the scalar DM and the Higgs, we impose a 
Z% symmetry (x — > — x) to ensure stability of the DM particle. The interactions in V(x, H) 
can also be constrained from WEP tests, and we will explore this in the next two sections 
for scalar DM. For vector-like fermionic DM no renormalizable couplings exist between the 
Higgs and DM. Such couplings can however be present for chiral DM. 

We consider the impact of a WEP violating force in the dark sector via the interactions 
of DM with the ultralight scalar <fi as in Eq.(12). These couplings are gauge invariant for 
fermionic DM only for vector-like gauge representations to which we restrict our attention. 
For chiral fermionic DM, the coupling to <fi can only arise from higher dimension operators by 
gauge invariance. Assuming no other low-energy degrees of freedom besides those of the SM 
plus the x and 0, the dark sector interactions (12) induce a coupling of to the SM fermions 
at two-loop order, as illustrated in Fig. 2. The graph involving only virtual U(l)y gauge 
bosons [left panel of Fig. 2] directly generate the operators Of in Eq. (48), while the one- 
particle irreducible diagram involving both W and B bosons in the 4>+f-^H + f "Compton 




(58) 
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FIG. 2: Two- loop diagrams in WIMP DM models that generate the operators Of in Eq. (48). 
Thus, after electroweak symmetry breaking the ultralight scalar couples to SM fermions. 

amplitude" [right panel of Fig. (2)] generates operators of the form Qi(p — Jp )Q, etc. As 
noted earlier, operators of this type can be expressed in terms of Of using the equations of 
motion, indicated symbolically by the presence of the H field on the external leg in the right 
panel of Fig. 2. In either case, the Wilson coefficients Cf are proportional to the Yukawa 
matrices due to the Higgs insertions. After the neutral component of the Higgs field obtains 
a vev, the loop-induced operators Of give rise to the interactions ff<f) of Eq. (16). Any 
mixing between the ultralight scalar and the Higgs will also contribute, corresponding to 
the second term as usual in Eq. (54). 

For SU(2)l triplet DM with hypercharge Y = 0, only the SU(2)i gauge boson exchange 
diagrams of the right panel of Fig. (2) contribute. The resulting coupling of <fi to the SM 
fermions are as in Eq. (16), with 

where we have employed naive dimensional analysis (NDA) to estimate the first term on the 
RHS of Eq. (59). Although the precise 0(1) coefficient C3 can be obtained from a complete 
computation, for our purposes of arriving at order-of-magnitude relationships between (3 
and 7], the NDA expression suffices. 3 We note that the sum of all loop graphs of the type 
in Fig. (2) is finite because we began with only renormalizable couplings and the operators 
Of have dimension n = 5. We also observe that the coupling to different species of fermions 
is universal since we have factored out the explicit dependence on the Yukawa coupling in 
the definition of the gj in Eq. (16). 

For SU (2) l multiplet DM with hypercharge F^O, the induced couplings of the ultralight 



3 The subscript in C3 refers to the dimension of the triplet representation of SU(2)^. 
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scalar to DM is given by 

2 , . 

1, . , 



gf ^ CN ^y^ A + CYY ^y^ A _ sin ^ _ (60) 



where Cn,y are 0(1) coefficients that, as before, can be obtained from a complete two-loop 
computation. We observe that the first terms on the RHS of Eqs. (59) and (60) are universal 
for different fermion species and come from the exchange of the SU (2) l gauge bosons W a 
in Fig. (2). The last terms containing sin 9 are also universal, having been generated from 
the mixing between the Higgs and light scalar <fi. The middle term in Eq. (60), involving the 
square of the SM hypercharge Y, are non-universal and are generated by the exchange of 
U(l) Y gauge boson B. We point out that such minimal WIMP DM models with non-zero 
hypercharge are typically ruled out [74] by direct detection experiments. Here we discuss 
these minimal DM models with non-zero hypercharge, only as illustrative examples keeping 
in mind that such DM could be part of a non-minimal extension which avoids the direct 
detection bounds. 

A similar analysis can be performed for other WIMP models of DM that may involve 
additional degrees of freedom. In supersymmetry, for example, the DM matter particle x 
is a linear superposition of winos, binos, and higgsinos. In addition there are squark and 
slepton particles which give interactions of the type Xipipx + h.c In theories with such a 
spectrum of particles one can induce a coupling of <fi to ordinary matter via virtual DM at 
one loop as shown in Fig. 3 4 . If the ultralight scalar <j) is the scalar component of a singlet 
superfield S, a superpotential term of the form (// + g x S)H u ■ H d will lead to a coupling to 
fermions of the form 

If x is primarily a bino, then A ~ gy, the hypercharge coupling. If x is primarily Higgsino, 
the coupling of (f> to the light quarks will be suppressed. The coupling of <p will be primarily 
to the top quark which has order one Yukawa couplings. Thus, in such models it is possible 
to induce a stronger WEP violating coupling to ordinary matter at one loop leading to 
bigger effects in Eotvos experiments. For the sake of brevity, we do not consider such non- 
minimal scenarios and we will only focus on minimal DM models without additional degrees 
of freedom such as squarks and sleptons. 

It is possible that the loop-induced Of operator contributions to Qf and those generated 
indirectly by H-S mixing (proportional to sin^) are individually much larger than gj yet 



4 Of course the presence of an ultralight scalar would introduce a new hierarchy problem which spoils the 
main motivation for supersymmetric theories. Here we invoke supersymmetry simply as a familiar example 
to illustrate the possibility of new types of interactions that can induce a coupling of the ultralight scalar 
to ordinary matter. 
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FIG. 3: DM-induced coupling of 4> to SM fermions at one loop in the presence of additional squark 
and slepton like degrees of freedom. 



cancel to produce a much smaller coupling. However, away from this special region of 
parameter space, each contribution will be roughly no larger in magnitude than gj itself (as 
already discussed following Eq.(54)). In this generic case, we are able to obtain expectations 
for the size of Eotvos parameters in our illustrative minimal WIMP dark sector models, for 
a given value of /3, from the two-loop gauge contribution. We will consider the contribution 
of H — S mixing to gj in the next section. We note that the contributions from the SU(2)i 
gauge bosons are generically an order of magnitude larger than those from the hypercharge 
gauge bosons due to the relative sizes of their couplings [leading to the additional factor 
of 1/16 in the second term of Eq. (60)]. Consequently, for purposes of making order-of- 
magnitude estimates, we may employ the expressions for t]s in the presence of universal 
couplings given in Eqs. (29) and (36) with 

2 - 
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■9 X £> 



(62) 



Expressing g x £ x in terms of j3 we then obtain 
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(63) 



(64) 



_(N p + N n ) + (m e /m N )N e 

Numerically the bounds in Eqs. (63) and (64) are (for Be and Ti laboratory samples with 
\Z 1 JA 1 - Z 2 /A 2 \ ~ 1/72) 

1dm > 1 A-10 



Ve 



> 10' 

> 10~ 16 



(65) 
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FIG. 4: An estimate of the allowed region in the (j] DME -,l3) parameter space for minimal WIMP 
DM models. The curves in the figures give an estimate of for rj DM _ for a given value of (3 from 
the two loop diagrams in Fig. 2. The shaded region is unlikely for typical WIMP models. Using 
the observational constraint (3 < 0.2, the allowed region is further restricted to the left of the 
vertical line. The estimates in the above figures for rj DME , for (3 < 0.2, are far below the current 
experimental bounds f] DM <■ 10 -5 , rj E ;$ 10 -13 . An improvement of about five orders of magnitude 
would be required in Eotvos experiments to fully probe the allowed parameter space for j3 = 0.2 
by measuring rj E . This is within reach of the MiniSTEP [59] proposal. 

which are shown in Fig. 4 as the allowed regions for typical miminal WIMP DM models in 
the (ti dme ,P) parameter space. The curve in the left and right figures gives as estimate of 
the minimum size for r\ DU and 7] E respectively as a function of (3. One can also estimate the 
ratio r) B /r) DM from Eqs. (63) and (64) to be approximately 

r) E /r] DM ~ 10- 6 . (66) 

For (3 = 0.2, marked by the vertical lines in Fig. 4, the current upper bound from 
galactic dynamics [54], we see that the lower bounds for typical WIMP DM models are 
i] DM > 4 x 10 -12 and r] E > 4 x 10 -18 . These lower bounds are far below the current 
experimental upper bounds shown in Eq.(10). An improvement of about five to seven orders 
of magnitude in Eotvos experiments would be required in order to probe these expectations 
of WIMP DM models. The MiniSTEP [59] experiment, which is currently under study, is 
expected to reach a sensitivity for r] E of about 10~ 18 and might be able to probe the lower 
bounds of these WIMP models. However, if (3 < 0.05 as indicated by a recent analysis [57] 
of the CMB and large scale structure formation, the lower bounds on rj DM E are far beyond 
current and future planned experiments. If an effect is detected in r] DME far above the 
expectations in Fig. (4) it would suggest the possibility that the coupling of (b to the SM 
fermions is mostly via h — <f> mixing corresponding to the last term appearing in Eqs. (59) 
and (60). One could extract a value for sin 9 and derive implications for various DM scenarios 
as discussed in the following section. Another possibility that might explain an effect above 
the expectation in Fig. 4 would be a stronger induced coupling of (b to ordinary matter 
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FIG. 5: Direct detection process for the scalar singlet and triplet x via t-channel Higgs exchange 
with a nucleon. The magnitude of the detection rate is determined by the strength of the coupling 
a,2- If X couples to the ultralight scalar cf>, the size of a 2 and thus the detection rate is bound by 
WEP violation constraints. 

in non-minimal DM models; for example a one loop coupling of <f) to ordinary matter (see 
Fig. 3) in the presence of additional squark degrees of freedom. 

VI. WEP TESTS, DIRECT DETECTION, AND HIGGS BOSON DECAYS 

As observed in Ref. [36], the presence of a non- vanishing (3 of astrophysically interesting 
magnitude, together with present limits on r\ B DM can imply upper bounds on the size of 
DM-nucleus cross sections relevant for direct detection experiments. Here we analyze these 
bounds in detail for the illustrative cases of scalar DM scenarios and argue that upper 
bounds on the DM-nucleus cross sections are less stringent than obtained in Ref. [36]. We 
further comment on the analysis of Ref. [36] at the end of section VII. We also consider the 
implications of a dark force for the DM relic density and derive corresponding constraints. 
Finally, using a light scalar triplet, as part of a multicomponent DM scenario, we show 
how the presence of a non- vanishing (3 - together with experimental limits on r] E DM - can 
preclude observable shifts in the rate for the Higgs boson to decay to two photons as one 
might otherwise expect. 

To include the full set of possible renormalizable interactions between the DM, SM fields, 
and ultralight scalar, we expand the scalar potential of Eq. (38), imposing the (x — x) 
symmetry need to prevent DM decays: 

V(H, S, X ) = V(H, S) + ^ Mix 2 + + a 2 H^H X 2 + 9 xX 2 S + X xsX 2 S 2 . (67) 

For the scalar singlet case, x is a rea l field, while for the real triplet with components x° 
and x ± y one nas [75] 

X 2 =(x°) 2 + 2xV ■ (68) 
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We take Mq and a 2 to be positive in order prevent a non-vanishing vev for \ an d the oc- 
currence of phenomenologically unacceptable cosmological domain walls. The experimental 
constraints on this DM model, for g x = 0, were recently explored in [60, 61]. 

After electroweak symmetry breaking, the H^Hx 2 term generates a contribution to the 
DM mass: 

Ml = Ml + a 2 v 2 . (69) 

Henceforth, we will take v = 246 GeV, M x , a 2 , and the mass of the SM-like Higgs boson 
(rrih) as independent parameters. All of them govern the x _nuc l eus cross section, whose 
leading order amplitude is generated by t-channel Higgs exchange as in Fig. 5 and is given 
by 



M^-^fNN , (70) 



m l 



where we have neglected the t-dependence of the amplitude for simplicity. Note that since 
the real triplet has zero hypercharge, the elastic DM-nucleus scattering has no contribution 
form Z-boson exchange at tree level. The corresponding cross section is 

gj qW^n (71) 

where, for simplicity, we have dropped the dependence on momentum transfer to the nucleus. 
[Recall that g h ~ 1.71 x 10~ 3 is the coupling of the Higgs to the nucleon as defined in Eq.(34)]. 
Note that the cross section decreases for increasing M x or decreasing a 2 . Note also that the 
coupling a 2 , together with the masses M x and m h , control the x Te ^ c density through the 
annihilation diagrams of Fig. 6. For M x ~ m h /2 for singlet DM or a light triplet in the 
multicomponent DM scenario, the Higgs exchange contribution becomes large, requiring a 
suppression of a 2 in order to maintain the observed CDM relic density. In what follows, we 
will generally avoid this regime. 



A. WEP Tests and Ultralight-Scalar-Higgs Mixing 



A relation between the Higgs-exchange contribution to ct x n and r] E DM arises for non- 
vanishing j3 because the parameters a 2 , M x , and mh that enter the cross section also control 
the strength of the DM-loop induced mixing between the Higgs boson and the ultralight 
scalar. After electroweak symmetry-breaking, these loop effects generate contributions to 
the mass-squared parameters \x\ is and The parameter /x| contributes only to the mass 
m| and n\ s contributes to sin# and m^. As with the contributions from Higgs loops to /x| 
discussed earlier in Section IV C, the DM loop contributions to this mass-squared parameter 
will also require the introduction of fine tuning to maintain a sufficiently small mass for the 
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FIG. 6: Annihilation diagrams for the scalar singlet DM. For the scalar triplet DM one has in 
addition the usual annihilation diagrams mediated by gauge interactions. If the gauge interactions 
of the triplet dominate the dynamics of annihilation, a DM mass of around 2 TeV is needed to 
saturate the relic density. 



ultralight scalar. Furthermore, as already mentioned and discussed in more detail in section 
VII, the finite renormalized parameter /i| must be further restricted in parameter space in 
order to maintain < 10~ 25 eV along with a dark force large enough to be observed. We 
implicitly assume that we are in this region of parameter space, conducive to the observation 
of a long range dark force. We will discuss the implications of other regions in parameter 
space in section VII. 

We begin by observing that in addition to the direct coupling 5± of S to the Higgs via 
the operator H^HS, a DM-induced 0-matter coupling arises from the one-loop contribution 
to this operator through the second diagram of Fig. 7. After renormalization in the MS 
scheme, the resulting finite coefficient of his operator is 

where 

^ = fl, singlet x, ^ 
1 3, triplet x, 

The factor of k = 3 appears in the case of the triplet x due to the three components 
of the triplet traversing the loop in the second diagram of Fig. 7. Here 5\(fi) is the finite, 
scale-dependent coupling counterterm whose numerical value is a priori unknown and whose 
presence is required to ensure renormalization group (RG) invariance of the physical prop- 
erties of the and h. Note that the mass parameter Mq (taken here to be positive) rather 
than M 2 appears in the argument of the logarithm since we are working in the theory before 
electroweak symmetry-breaking. 

We also observe that the x 2 S interaction will generate a contribution to the mass param- 
eter /i| as it yields the non-vanishing contribution to the S self energy: 



2 



J - 7 + In 4vr + In /i 2 - F (M 2 , M 2 , p 2 



(74) 
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FIG. 7: One loop diagrams which which contribute to the effective potential V(H, S) with one 
external S field. After electroweak symmetry breaking the effective potential contributes to Higgs- 
ultralight-scalar mixing. The first two diagrams are UV divergent and contribute to the renormal- 
ization of the S-tadpole and the coupling 5\ respectively. The remaining diagrams mix into higher 
dimensional operators and give a finite contribution to Higgs-ultralight-scalar mixing as explained 
in appendix B. 

As with the case of the logarithmically divergent Higgs contribution T>(p 2 ) H i HS of Eq. (55), 
the DM-loop contribution to the self energy requires a corresponding /i-dependence in Sfj, 2 s (fi) 
to maintain RG invariance of the pole mass that governs the range of the dark force. Large 
DM-loop contributions to /i| can be minimized at all scales by taking g x to be sufficiently 
small: g x < 4Trm s ° le . Doing so, however, would preclude a value of (5 of astrophysically 
relevant strength. Alternatively, one may allow for a much larger, phenomenologically in- 
teresting magntiude for g x and maintain a small /i| by invoking fine-tuning between the 
one-loop contribution of Eq. (74) and Sm 2 (//). 

A similar set of alternatives applies to the renormalized coupling S[ en . One could require 
that the product g x a<i be sufficiently small in magnitude, with a correspondingly tiny Si(fi), 
so that the induced H-S mixing is consistent with the present bounds on t/e,dm- To obtain 
a large (3, one must then take a 2 to be sufficiently small, implying an upper bound on the 
Higgs exchange contribution to the DM-nucleus cross section via Eq. (71). This choice is 
essentially the strategy followed in Ref. [36] to obtain upper bounds on u x n- However, as 
seen in Eq.(41), Si contributes to fihs and thus to the mass via Eq.(47). The condition of 
m<£ < 10~ 25 eV gives a much stronger naturalness constraint on a<i forcing it to be essentially 
zero for a non-zero dark force. The constraints from i] E DM are thus not relevant in such a 
naturalness analysis. We will also discuss this in more detail in section VII. 

In what follows, we will instead allow for fine tuning in both Si since we have already 
allowed fine tuning for [ig. We show that assuming /i| is restricted in parameter space to 
satisfy < 10~ 25 eV for any value of fi\ s in Eq.(47), we can obtain upper bounds on a X N 
by analyzing finite, one-loop contributions to fi\ s , from higher dimensional operators after 
EWSB, and their implications for WEP tests. The other regions in parameter space and 
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their implications will be discussed in section VII. To that end, consider the third diagram 
of Fig. 7, which generates a contribution to the dimension five operator 

C 2 (H^H)(H^H)S, (75) 

where in the coefficient C 2 is finite and given by 

Since this contribution is finite there is no counterterm involved in determing the value 
of C 2 . After electroweak symmetry breaking this term will generate a contribution to the 
off-diagonal elements in the h-S mass-squared matrix 

fi 2 hS = 2C 2 v 3 + 5 lV , (77) 

leading to an h-S mixing angle 6 

tang - _ vis _ 2C 2 v* + Sjv 

which was defined in Eqs. (42) and (45). Since this mixing implies a coupling of <fi ~ S 
to matter, the loop-induced coefficient C 2 will contribute to the Eotvos parameters rj DME . 
Given the dependence of C 2 on a 2 and the absence of any fine-tuning in this parameter, we 
obtain an upper bound on <j x n for non-vanishing f3 as described below. 

Before doing so, we observe the contribution to [i 2 hS from full series of diagrams appearing 
in Fig. 7 (plus the tadpole graph generated by the \ 2 S interaction) can be evaluated in a 
straightforward way as outlined in the appendix B. After renormalization, the result is 



M. 



An 2 \ fj, 



,2 



,3 



k 9 ^ — (79) 



Apart from an overall constant in the first term and the replacement M — > M x , this 
expression is the same as we obtained using the contributions to the H^HS and (WH) 2 S 
operators from the second and third diagrams of Fig. 7. The expression in Eq. (79) has 
the advantage that it depends on the tree-level x mass after electroweak symmetry breaking 
rather than on the parameter Mq as in the effective operator analysis. We will henceforth use 
the finite, second term in Eq. (79) to derive an upper bound on Higgs exchange contributions 

to (J x n ■ 

To that end, we write the mixing angle as 

(I 2 Q S Yen v a 2 S Yen v 

4tt 2 M 2 m 2 m 2 vr 3 / 2 M P m 2 m 2 v ' 

where 5l en denotes the quantity in square brackets in Eq. (79). The mixing angle sin# 
also characterizes the universal H-S mixing contribution to the Eotvos parameters f] EDM . 
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FIG. 8: Upper bound on 02 in the singlet (red-solid) and real triplet (black-dotted) scalar DM 
models as a function of rj DM /(3 2 (left panel) and r] E /f3 2 (right panel). We have used = 120 
GeV and assume (3 = 0.2 to discuss the resulting bounds on 02 from Eotvds experiments. In 
the top left and right plots, the vertical black lines on the right correspond to the upper bounds 
Vdm < 10 -5 and rj E < 10 -13 respectively. These vertical black lines will move to the left with 
further improvements in Eotvds experiments as indicated by the left-pointing arrow in each plot. 
The bottom left and right plots show the region closer to the expected future bounds, from the 
MiniSTEP experiment, of r\ DM < 10 -10 and r\ B < 10~ 18 respectively as indicated by the vertical 
black line in each plot. We explore the implications of these bounds on 02 from Eotvds experiments 
for specific observables in sections VIC and VI D. 



We now require that the contribution from each term on the RHS of Eq. (80) to r] E DM be 
no larger than the experimental limits on these parameters. As discussed previously, the 
different parametric dependence of each term and avoiding slices of parameter space with 
unnatural cancellations between the two allows us to treat each one separately. Considering 
only the first term proportional to a|, using Eqs. (29), and (36) with 

• a mN 

g — > — sinf (81) 
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and expressing g x in terms of (3, we obtain 




' g h (N p + N n ) + (m e /v)N e ~ 
X _(N p + N n ) + (m e /m N )N e _ 

Eqs. (82) and (83) imply that for fixed (5 and m/,, the experimental bounds on r] E DM 
translate into bounds on a 2 as shown in Fig. 8. The solid red curves and the dashed black 
curves correspond to the bounds on a 2 in the singlet x an d real triplet x models respectively. 
For f3 = 0.2 the current bounds in Eotvos experiments imply t] DM /(3 2 < 2.5 x 10 -4 and 
f] E /P 2 < 2.5 x 10~ 12 respectively. The tighter bound from r\ E implies a 2 < 0.042 and 
a 2 < 0.026 for the singlet and triplet x models respectively. The possible future experiment 
like MiniSTEP is expected to improve the the sensitivity oir) DM E by five orders of magnitude. 
From Eq.(83) this would tighten the bound on a 2 by an additional factor of lO 5 / 4 for a non- 
zero (5. 

B. WEP Tests and Scalar Singlet DM Relic Density 

For the scalar singlet DM, the DM relic density is determined entirely by the parameter 
a 2 for fixed DM and Higgs masses. This feature can be seen from Fig. 6, where a 2 enters the 
amplitude for each annihilation diagram and, thus, determines the DM annihilation rate. In 
particular, the value of a 2 must be sufficiently large so that the DM relic density does not 
over saturate the observed value. Thus, the requirement that the singlet DM relic density 
smaller than the total DM relic density, 0, DM < &dm, leads to a minimum value a™ m for 
fixed values of M x and m h . If the WEP bounds on a 2 imply that a 2 < a r 2 m , then a dark 
force of the corresponding strength will be ruled out. 

In order to illustrate this interplay, we refer to Fig. 3 of [61]. The parameters a 2 and M x 
are the same as A and mo respectively in the notation of [61]. From Fig. 3 of that work, 
we see that for DM masses in the 0-50 GeV range, the required value of a 2 ranges from 
~ 0.16 — 0.05 respectively for = 120 GeV. On the other hand, from the upper curve in 
the top right graph of Fig. 8, we see that a 2 ^ 0.045 for (3 = 0.2 from the current bound 
of rj E < 10~ 13 (vertical black line at right). This WEP constraint a 2 < 0.045 implies an 
over-density of DM in the range < M x < 50 GeV thus ruling out the possibility of a dark 
force with f3 > 0.2 in this range of parameter space. We give sample points in the parameter 
space of singlet DM models in Table II. 

A more detailed analysis can be performed to rule out even smaller values of (5 depending 
on the DM mass in the scalar singlet model. Future Eotvos experiments with the sensitivity 
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G2relic 


M x (GeV) 


Expectation for^f- 


(3 = 0.2 


0.15 


20 


4 x 10- 10 


Excluded 


0.10 


40 


7 x 10~ n 


Excluded 


0.02 


100 


1 x 10~ 13 


Allowed 



TABLE II: The first two columns give sample points in the (a2 re ii c , M x ) space of scalar singlet DM 
models with a Higgs mass of m/, = 120 GeV. The third column gives an expectation for r] E /(3 2 
from Eq.(83). The fourth column uses the current bound of rj E < 10~ 13 to determine whether a 
dark force of (3 = 0.2 is ruled out. One can equivalently compare the different values of cyclic with 
the WEP bound on 02, at (3 = 0.2, in top right graph of Fig. 8 at the far right verticle line. 

of MiniSTEP [59] which are expected to reach a sensitivity of r\ E < 1CT 18 , could require 
bound of a 2 ^ 0.0025 for (3 = 0.2 as seen in the bottom right graph of Fig. 8. In this case, 
one can rule out (3 < 0.2 even for DM masses above 60 GeV which require smaller values 
of a 2 in order to get the right relic density. As seen in Fig. 3 of [61], larger values of the 
Higgs mass typically imply much larger values of a 2 . For example, a Higgs mass of 200 GeV 
implies a range of a 2 of ~ 0.42 — 0.05 for the DM mass range of — 80 GeV thus ruling out 
the possibility of (3 > 0.2 in order to prevent an over-density of DM. Thus, the bound on a 2 
from WEP constraints is a powerful probe of a dark force in the scalar singlet DM model. 

For the scalar real triplet DM model, the DM relic density is determined by gauge inter- 
actions in addition to the parameter a 2 . In this case, the WEP bound on a 2 shown in Fig. 8 
does not necessarily rule out a dark force since the correct relic density can still be obtained 
from annihilation diagrams that proceed via gauge interactions that are independent of a 2 . 
For example, the bound of a 2 < 0.02 implied by r] E < 10~ 13 for (3 = 0.2, as shown in the 
top right graph of Fig. 8, implies that the annihilation rate will be dominated by gauge 
interactions. 



C. WEP Tests and DM-Nucleus Cross-Sections 

The current bounds on a 2 for a non-zero (3 in the dark sector, will also lead to upper 
bounds on the Higgs exchange contributions to the direct detection cross-section. From Eq. 
(71), the parameter a 2 can be written in terms of the tree level cross-section ct x n, which 
proceeds via a t-channel higgs exchange, as 

a\ = 



ir(M x + mjv) 2 m^ 
2 9 2 

9h V m N 



a 



(84) 
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Substituting into Eqs. (82) and (83), denning the quantities 



F = 



E 



,Z 1 ^ 2 \f 



m N 



m„ 



x, 



m N 



>,p x q,n)} 



g h {N p + N n ) + {m e /v)N e 
(N p + N n ) + (m e /m N )N e 



9h, 



(85) 
(86) 



we obtain the following relations between (M x + rriN) 2 a X N and the Higgs-exchange contri- 
butions to the DM-nucleus cross section: 



(M x + m N ) 2 a xN 



(M x + m N f a xN 



Higgs exch 



Higgs exch 



18 
7^ 

6"7T 



9h 



1 V M 



F (3 2 



9h 



( m N \ 2 






V2 1/2 
<E 


\m h J 


[(v) 


FE_ 


P 



(87) 
(88) 



The experimental limits on the Eotvos parameters, together with the foregoing expres- 
sions, lead to bounds on the Higgs exchange contributions to the DM-nucleus cross sections. 
These bounds can be brought into the numerically convenient form as a function of rj DM //3 2 
for the singlet x as 



M x + m N 
100 GeV 

and for the triplet x as 



S < ( L1 >< 1Q4 )^ 



"100 GeV" 
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(89) 
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(90) 



Similarly, the bounds as a function of rj E /j3 2 can be brought into the numerically convenient 
form for the singlet x as 



M x + rriN 
100 GeV 



a xN 
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< 8.4 x 10 2 g 2 h 
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(91) 



and for the real triplet x as 



M x + rriN 
2TeV 



°XN_ 

lpb 



< 0.7 gl 



100 GeV 



Zi 
A 1 



Zi 
A 2 



-1/2 „l/2 

'E 

(3 ' 



(92) 



where 
upper 



Z2 

A 2 



1/72 for Beryllium and Titanium samples in Eotvos experiments. The 
Dounds for the Higgs exchange tcontribution to the direct detection cross section of 
the singlet (left panel) and real triplet (right panel) DM in the presence of a dark force, as 
determined by Eqs. (89), (90), (91), and (92), are shown in Figs. 9 and 10 as a function of 
t] dm //3 2 and rj E /j3 2 respectively. In each graph we show three sample curves corresponding to 
Higgs mass choices of m h = 120, 130, and 140 GeV as indicated. In the top row of Figs. 9 and 
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FIG. 9: Upper bounds on the Higgs exchange tree level direct detection cross section of scalar 
singlet (left panel) and real triplet (right panel) DM implied by a dark force as a function of rj DM / 1 1 '■ 
To be specific, we assume (3 = 0.2 and discuss the implied bounds. In the top left and top right 
plots the vertical black line on the right corresponds to rj DM < 1CP 5 . In the bottom left and bottom 
right plots, the vertical black line corresponds to the expected future sensitivity of r] DM < 10~ 10 . In 
all plots, the three lines from to bottom correspond to the bounds for the Higgs masses of 120, 130, 
and 140 GeV respectively. The size of these bounds compared to current and future sensitivities 
for direct detection experiments is discussed in the text. 

10, the vertical black lines on the right correspond to the current WEP bounds of r] DM < 10 -5 
and r] E < 10~ 13 for our benchmark value of (3 — 0.2 which we use throughout this discussion. 
In the bottom row of Figs. 9 and 10, the vertical black lines labeled 'MiniSTEP' correspond 
to the expected sensitivity of 7/ DM ~ 10~ 10 and r] E ~ 10~ 18 from a possible future experiment 
like MiniSTEP [59]. We see that the bounds on the Higgs exchange contribution to the DM- 
nucleus cross-sections are typically much stronger from WEP violation constraints on r] E 
compared to those on r\ DM . However, since the DM-nucleus cross-section bounds depend 
linearly on r\ DU and on the square root of rj E , with enough improvement the bound from 
r\ DM could become stronger. In the following discussion, we focus only on the direct detection 
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FIG. 10: Upper bounds on the Higgs exchange tree level direct detection cross-section of scalar 
singlet (left panel) and real triplet (right panel) DM implied by a dark force as a function of rj E /P 2 . 
To be specific, we assume (3 = 0.2 and discuss the implied bounds. In the top left and top right 
plots, the vertical black lines on the right correspond to the current bound of r\ E < 10 -13 . The 
the bottom left and bottom right plots, the vertical black line corresponds to the expected future 
sensitivity of r\ E = 10~ 18 . In all plots, the three curves from top to bottom, correspond to the 
Higgs masses of 120, 130, and 140 GeV respectively. The size of these bounds with current and 
future sensitivities of direct detection experiments is discussed in the text. 

bounds from rj E shown in Fig. 10. 

For the scalar singlet DM, the DM-nucleus cross-section bound from WEP tests does not 
yield any more information than the bound on which has already been discussed. This 
is due to the fact that a-i determines the DM matter relic density entirely for fixed DM and 
Higgs masses. If the WEP violation bound on 02 is too strong, the resulting DM relic density 
will be too large over-closing the universe and thus ruling out the dark force. The bound 
on the DM-nucleus cross-section resulting from the corresponding WEP violation bound on 
a 2 , is thus not useful since it is already ruled out. 

However, the bound on the scalar singlet DM-nucleus cross-section can be useful in con- 
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M x = 50 GeV 


Experiment 


Sensitivity 
a xN (pb) 


Sensitivity 

r 1 2 r i 
I 100 GcV J [ 1 pb 




CDMS [76] 
XENON10 [17] 
CDMS (2007 [77]) 
WARP (140 kg) [78] 
SuperCDMS (Phase A) [79] 
WARP (1 ton) [80] 


1.6 x 10~ 7 
4.5 x 10~ 8 
1 x 10~ 8 
3 x 10~ 8 

1 x 10~ 9 

2 x 10- 10 


4.1 X 10~ 8 

1.2 x 10~ 8 
3 x 10" 9 
8 x 10" 9 
3 x 10~ 10 
5 x 10- 10 



TABLE III: Sensitivities for DM direct detection cross sections in different experiments. These 
sensitivities are for 50 GeV DM corresponding to the most sensitive mass window. We see that the 
XENON10, CDMS (2007), WARP (140 kg), SuperCDMS, and WARP (1 ton) experiments have 
enough sensitivity to probe the bounds on the direct detection cross sections in Fig. 10 for singlet 
DM coupled to a WEP violating force. 

straining the size of j3 in a multicomponent DM scenario where the scalar singlet is only a 
fraction of the DM. For larger values of /3, as already discussed, the bounds on a 2 from WEP 
tests are too strong leading to an over-closed universe. For smaller values of (3 the bound 
on a 2 becomes weaker as seen from Eqs. (82) and (83). For small enough values of /3, the 
upper bound on a 2 would be consistent with an under-relic-density of the singlet scalar. A 
multicomponent DM scenario can also have a 2 consistent with an under-relic-density for the 
scalar singlet and in this case the WEP violation constraints on a 2 can lead to interesting 
bounds on the DM-nucleus cross-section. 

For the scalar real triplet DM, the DM relic density is determined by a 2 and gauge 
interactions in general. However, the tree level DM-nucleus cross-section proceeds only 
via a t-channel Higgs exchange and its size is determined by a 2 . We point out that the 
bound in Fig. 10 constrains the tree level Higgs exchange diagram but not the one loop 
diagrams, which proceed via gauge interactions and the Higgs coupling to the nucleus and 
is independent of a 2 . Thus, if the observed DM-nucleus cross-section is of the size explained 
by this one loop diagram a dark force cannot be ruled out. 

In Table III we show the sensitivities of current and future DM detection experiments, 
taken from Table I of [60]. We see from Fig. 10 that it will be difficult for current and future 
direct detection experiments to probe the upper bound on the DM-nucleus-Higgs-exchange 
cross-sections, for scalar triplet DM, for values of (3 that are astrophysically interesting 
allowing one to rule out this possibility. One would need a significant deviation from the 
expected cross-section from one loop gauge diagram, indicating a large value of a 2 , to rule 
out a significant dark force. For smaller enough values of /3, the DM-nucleus cross-section 
bounds should be within reach of current or future experiments. The bounds we have 
derived on the DM-nucleus cross-sections are much weaker than those in [36] since our 
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FIG. 11: Contributions to the h — > 77 rate from virtual x^ loops. 

analysis constrains higher dimension operators with finite coefficients while the work of [36] 
had to rely on naturalness arguments to constrain renormalizable couplings. As we have 
discussed earlier, since the ultralight scalar mass is itself fine tuned we have avoided using 
naturalness arguments. 

D. WEP Tests and Higgs Decays 

We have shown in the last section that WEP constraints lead to upper bounds on the tree 
level DM-nucleus cross-sections for the scalar singlet and real triplet x models. However, 
if the dark sector is made up of a rich spectrum of DM particles of different species, direct 
detection of any species that makes up only a tiny fraction of the relic density becomes 
difficult. One example of such a DM species is the neutral component of the real triplet 
scalar \ with a mass far below a TeV. For masses below 500 GeV, the triplet DM will 
make up less than 10% of the relic density [74]. The astrophysical effects of a dark force 
experienced by such a species would be too small to be detected. In this section, we show 
that when direct detection experiments or astrophysical observations fail to constrain dark 
forces, collider signals might still harbor information on dark forces. Fig. 2 of [75] shows 
the size of the shift in the h — > 77 rate for typical values of the parameter 02- We have 
reproduced this figure as shown on the left in Fig. 12. We plot the quantity 

<H % ) = 100 x T SM {h - 77) > ( 93 ) 

For specificity we focus on the real scalar triplet x discussed in the last section, but 
with a mass less than 200 GeV, and examine the implications of a dark force on collider 
signals. The analysis of [75] showed that one potential signature of the scalar triplet would 
be a modification of the h — > 77 decay rate due to the virtual charged components of 
the x triplet traversing the loop shown in Fig. 11. In the rest of this section we focus on 
this channel. For a heavy Higgs, a similar analysis can be done for h — > jZ, ZZ, W + W~ . 
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FIG. 12: The left plot shows the deviation of the h — ► 77 rate compared to the SM prediction for 
typical values of the parameter a 2 as a function of the triplet mass M x . The right plot shows the 
magnitude of the allowed shift in the h — > 77 rate in the presence of a dark force. The bound on 
this allowed shift arises due to the bound on a 2 from WEP violation constraints on t/e as seen in 
the top right plot of Fig. 8. The typical values of a 2 in the left plot above, which lead to sizable 
deviations in the h — * 77 rate, are too big to be compatible with WEP violation constraints. 

As already discussed, WEP constraints imply an upper bound on the parameter 02 which 
determines the size of the contribution of Fig. 11 to h — > 77. The WEP bound on a 2 
translates into a bound on 5(%) which is shown in the right panel of Fig. 12 for different 
values of M x . Comparing the left plot of Fig. 12 with the top right graph in Fig. 8, we see 
that the current bounds on a 2 from r] E for a non-zero (3 can give non-trivial bounds on 5(%) 
that can be tested in colliders. The right plot in Fig. 12 gives the upper bound on |5(%)| as 
a function of tje/ ' (3 2 - For (3 = 0.2, we have the bound i] E /j3 2 < 2.5 x 10 12 coming from the 
current bound of rj E < 10~ 13 . We see that the bound on 5(%) for a dark force of f3 — 0.2 is 
well below one percent. Thus, any observed shift in h — > 77, that cannot be explained by 
physics observed at colliders and unrelated to Xi requires a significant contribution from the 
X loop implying a value for {3 much smaller than 0.2. If a non-zero value of a 2 is extracted 
from a study of h — > 77 decays, one can estimate the size of Vdme/0 2 from Eqs. (82) 
and (83) respectively and use the current bounds on r] DME to constrain the size of (3. For 
example, using = 120 GeV and the current bound of r] E < 10~ 13 , the non-zero values of 
a 2 = y^F, 1.0, 0.5 would imply that (3 < 7 x 10~ 5 , 2 x 10~ 4 , 9 x 10" 4 respectively. 
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VII. DARK FORCE PARAMETER SPACE 

Before concluding, we generally discuss the regions in parameter space of SM + \ + 
type models that are likely to give rise to an observable dark force. In particular, we discuss 
how the requirement < 10 _25 eV, neccessary to allow a dark force of intergalactic range, 
restricts the allowed parameter space. Recall that after EWSB and diagonalizing the mass 
matrix, the ultralight scalar mass is given by 

We showed in section VIA that fi 2 hS , which determines the mixing between the ultralight 
scalar and the Higgs, receives finite contributions from higher dimension operators whose 
size we constrained from WEP tests. In Eq.(80), the second term is the contribution to the 
mixing angle from the operator WHS after EWSB and the first term is the finite contri- 
bution from the sum of higher dimension operators induced via DM loops (see Fig. 7). The 
parameter /x| similarly receives finite contributions from higher dimension operators. For 
example, attaching one extra external S field to the DM loops in Fig. 7 will generate a tower 
of operators that contribute to fi s after EWSB. The lowest dimension non-renormalizable 
operator that contributes to /i s will be 

D 2 WHH^HS 2 , (95) 

with finite coefficient D 2 which can be estimated from NDA as 

D > ~ Jr 32 (96) 

We could sum the contribution of the entire tower of operators to pig as we did for the case 
of n 2 hS . However, the explicit sum is not needed for the following discussion. 

The requirement that < 1CT 25 eV now imply three types of possible regions in pa- 
rameter space: I, II, and III. We discuss each of these regions in turn below and relate them 
to the analysis of previous sections. 

Region I: In the first region of parameter space, there are no intricate cancellations of 
any kind among the terms in Eq.(94). Each term that goes into determining /i| and ix\ s is 
required to be of the order of m 2 . In this case, we can obtain an an approximate bound on 
a 2 from the D 2 which contributes to /i| ~ m| as 

4tt M 2 p m\ 3 x 1(T 39 
]3 2 ^u 2 ~^u 2 ~ < p 2 

For any observable non-zero value of (3, the above bound essentially forces a 2 to be zero. As 
already discussed, such a small value of a 2 will lead to an over-relic-density of scalar singlet 
DM over-closing the universe and is ruled out. Similar arguments can be made for WIMP 
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DM in which case D 2 will receive contributions that depend only on the gauge couplings 
and /3, thus ruling out any observable value for (3. In short, this region of parameter space 
is incompatible with the existence of a phenomenologically interesting dark force and thus 
not considered in this paper. 

Region II: In this region, /i| of Eq.(94) is chosen such that the condition < 1CT 25 
eV is always satisfied regardless of the size of n\ s and any finite contributions to it from 
higher dimension operators after EWSB. Furthermore, in this region there are no intricate 
cancellations between terms that determine fj,\ s so that one can put bounds on these from 
WEP tests. This is phenomenologically the most interesting region and was the focus of 
this paper. 

Region III: Finally, the third region corresponds to the case where there are intricate 
cancellations among various terms in /i| and n\ s individually. If there are intricate can- 
cellations between large terms in n 2 hS) then we cannot extract meaningful bounds on a 2 or 
the Higgs exchange direct detection cross-section from WEP tests. In this special region 
of parameter space the bounds on a 2 and the Higgs exchange contributions to the direct 
detection cross-sections derived in this paper do not apply. 

The analysis of Ref. [36] assumed no fine tuning of 5\ or n\ s against radiative corrections 
sensitive to the cutoff. As seen in Eq.(94), contributes to the light scalar mass and in 
the absence of fine tuning the strongest bound on a 2 comes from < 1CT 25 eV forcing a 2 
to be essentially zero. The bounds from WEP tests are relatively far weaker and not relevant 
for a 2 or equivalently for the DM-nucleus cross-section via Higgs exchange. Our analysis 
differs in that we allow for fine tunning in all renormalizable parameters, and then examine 
the different regions in the space of these renormalized parameters and the corresponding 
implications for terrestrial experiments. 

VIII. CONCLUSION 

The existence of a new long-range WEP-violating attractive force in the dark sector, a 
dark force comparable to gravity, can have interesting cosmological consequences, including 
an accelerated rate of structure formation and an explanation of certain features of the 
DM distribution and other astronomical observations [44, 48, 53]. Strong constraints for 
such a dark force comes from a study [54, 55] of the dynamics of satellite galaxies and 
the evolution of density perturbations constrained [57] by the CMB spectrum. The current 
bounds indicate the strength of a possible dark force to be less than < 20% of gravity from 
galactic dynamics and less than 5% of gravity from the CMB spectrum. 

We considered the consequences of such a dark force for terrestrial experiments. Ordinary 
matter will feel the effects of a dark force via virtual DM as long as the DM candidate is 
not sterile. Depending on the DM model, a dark force can lead to constraints on Eotvos 
experiments, DM-direct-detection experiments, and Higgs decay properties to be studied 
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at future colliders. We studied several minimal models of DM to illustrate the prospective 
implications of an astrophysically relevant dark force for terrestrial experiments. 

• We derived lower bounds on the size of the Eotvos parameters rj DM E for a non-zero dark 
force for minimal DM models. These Eotvos parameters measures the effect of WEP 
violation on ordinary matter which can arise through virtual DM that communicates 
the dark force to ordinary matter. We find that for light scalar singlet DM, relic 
density considerations and the experimental limits on r\ E rule out a dark force having 
strength of 20% of gravity, in large regions of parameter space. Future experiments 
with improved sensitivity could probe a dark force of this magnitude for heavier singlet 
DM. For minimal WIMP DM, the expected magnitudes of rj DME lie well below current 
and prospective sensitivities of terrestrial Eotvos experiments, but could be probed in 
a satellite-based experiment having the sensitivity of the MiniSTEP proposal. In 
non-mimimal WIMP DM models, it is possible to generate larger effects for Eotvos 
experiments that could be detected by the Microscope experiment. 

• WEP tests imply constraints on Higgs-exchange contributions to the DM-nucleus 
cross-sections. For scalar singlet DM, these bounds apply to the entire cross section. If 
the scalar singlet DM saturates the DM relic density, these bounds on the DM-nucleus 
cross-section do not give any information beyond the implications of WEP bounds 
on the DM relic density. If the scalar singlet does not saturate the DM relic density, 
allowing for stronger interactions with the Higgs, WEP tests provide useful bounds on 
the DM-nucleus cross-section. For WIMP DM, the WEP constraints on the Higgs-DM 
interactions give upper bounds on the contribution of DM-nucleus scattering via Higgs 
exchange, to the total cross-section. 

• For singlet DM, the current WEP bounds on the DM-nucleus cross-section in the 
presence of a dark force that is 20% of gravity, are typically within reach of current 
and future direct detection experiments. For scalar WIMP DM with a mass in the 
TeV range, the corresponding bounds on the DM-nucleus cross-section are typically 
beyond the reach of current and future direct detection experiments. If these scalar 
WIMPs are detected, it will rule out a dark force greater than 20% of gravity, implying 
a tighter upper bound on the dark force. 

• The WEP violation constraints on the DM-Higgs interactions can lead to constraints 
for collider physics. As a specific example, we derived testable bounds on the allowed 
shift in the h — > 77 rate when the Higgs couples to WIMP DM in the real triplet 
scalar representation with a mass less than 200 GeV. Such a light triplet will only con- 
tribute a tiny fraction of the DM relic density and could be part of a multicomponent 
DM scenario. The implied bounds on the Higgs to two photon rate or a dark force 
comparable to gravity, is far below the sensitivity of the LHC or the ILC. An observed 
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shift in h — > 77 attributed to the charged components of the triplet would rule out 
the dark force. 

• An observable scalar dark force with intergalactic range implies restrictions in the space 
of the renormalized parameters of the theory. These parameter space restrictions apply 
after the usual fine tuning of parameters against radiative corrections sensitive to the 
cutoff. 

Apart from these experimental implications, a notable theoretical consequence of an astro- 
physically interesting dark force mediated by an ultralight scalar is the need for substantial 
fine-tuning to preserve its tiny mass (m^ < 10~ 25 eV). In the DM scenarios considered here, 
divergent loop contributions associated with the DM or SM particles that interact with <fi 
would generate large contributions to that must be removed by fine-tuning unless the 
strength of the dark force is imperceptibly small. The discovery of such a dark force would 
introduce yet another mass hierarchy problem in particle physics. In our analysis, we have 
taken this need for fine-tuning at face value and have attempted to apply it consistently to 
the derivation of implications for terrestrial experiments. These consequences imply that 
direct detection experiments - together with Eotvos experiments and astrophysical observa- 
tions of satellite galaxies and structure formation - can be employed as part of a multifaceted 
probe of a long range force in the dark sector. 
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APPENDIX A: ATOMIC CHARGE TO MASS RATIO UNDER A COUPLING TO 

ULTRALIGHT SCALARS 



Here we give details for obtaining the atomic charge to mass ratio for a coupling to 
ultralight scalars. As discussed in the text and shown in Eq.(20), we can obtain the atomic 
coupling to <fi from the couplings of to SM particles by a matching calculation. From 
Eqs.(19) and (20), the coupling atomic coupling g A is given by 



Ug A = —{A\ ( 5> s m ff 9<z + J>m^ ) \A) + c g (A\G^G%\A) + c,(A\F^F, u \A)- 

(Al) 



m p 
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Note that with the normalization factor £4 the above equation is dimensionally consistent, 
since for fermionic atoms qa is dimensionless and for scalar atoms it has dimension one. 
Next we exploit the properties of the energy-momentum tensor in the low energy effective 
theory. The trace of the energy momentum tensor is given by 

K = Ir^cr + § v + E m w + E m &> ( A2 ) 

A93 AC q i 

where in this convention the QCD and QED beta functions at one loop are 

11 n (n\ n fn (at ^l dl a _ e 3 



3 C 2 (G) - -fC 2 (N c 



The mass of the atom A is given by 
Ma = (A\6#A) 



16tt 2 ' ye 12n 2 



Pe = ^~ 2 n f . (A3) 



(^Tr[G„GH + § \A) + (A| (j>^ + X>^) l A >- 



(A4) 

The atomic mass can also be expressed as 

Ma = Zrrip + (A — Z)m n + Zm e — £a, (A5) 

where Sa is the binding energy of the atom A. From here one can write the derivative with 
respect to the quark and electron masses as 

dM A „ dm p dm n d£ 

m q - = Zm q - \- (A - Z)m q - m„- — , 

dm q dm q dm q dm q 

m f - — = Zm e d e p - mp- — . (A6) 
dmp dmp 

We now use Eq. (A4) for M A and the Feynman-Hellman theorem to calculate f^q/j^^ and 
equate with Eq. (A6), we obtain 



dM A ... _ . . dm p dm n d£ A 

= {A\m q qq\A) = Zm q - \- (A - Z)m q - m q - — , 

dm n ■ ( i m ( j bm ( j bm 



in,,— = {A\m q qq\A) = Zffi q — - + (A - Z)m„- — - /// 

^q 

m ^^A = ^ me gg\A\ = Zm e 5 ei - m^. (A7) 
dmp dmp 

Using the relations of Eq.(A7) in Eq. (A4) we obtain an expression for the atomic matrix 
element of the gluon operator as 

/?3 r ,, ., / \ - dm p \ dm 



(A\^Tr[G„ u Gn\A) = M A -z(m e + J2^ q ^)-(A-Z)J2 

,— -(A\^ 
dm k 2e 



+ E-^-^i§^^)- ( A8 ) 
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Using Eqs.(A7) and (A8) in Eq. (Al) we finally arrive at the expression for £a9a'- 



U9A = M A + 

P3 



n k 



+ k^IF^IA), (A9) 

where we have introduced the index k which runs over the light quarks q and the charged 
leptons £ and the parameters (k and k are given by 

Ck = ^~~ IT "' K- = Cj 9 j^c g . (A10) 

We now utilize the expression for the nucleon mass in terms of the nucleon matrix element 
of the trace of the three flavor QCD energy momentum tensor 

m N = (N\0*\N) = (N\ ^-G^G? + J2m q qq \N), (All) 

the variation of the nucleon mass with respect to the mass of of a quark of flavor q is given 
by 

= (N\m g qq\N). (A12) 

Once again we have used the non-relativistic normalization of nucleon states. The nucleon 
matrix elements on the RHS are extracted from pion-nucleon scattering data using chiral 
perturbation theory. Experimentally their values are determined to be [68, 69] 

x u ,p = = (p\uu\p) ~ 0.019 m p /m u , 



x d, P = ^T d = (PIMP) ~ °- 041 m p /m d , 
x *,p = £ = (p\^ s \p) ~ 0.14m p /m s , 
x u,n = |fe = (n\uu\n) ~ 0.023 m n /m u , 
= ^z = (n\dd\n) ~ 0.034 m„/m d , 



dm. 



(n\ss\n) ~ 0.14 m n /m s . (A13) 



These numbers are taken from table 6 of [69]. In general rrikdSA/ ' dm^ and (AjF^F^IA) are 
not analytically calculable, at least for large atoms, and will contribute to the uncertainty 
in the atomic charge to mass ratio. The atomic charge to mass ratio can be finally written 

as 



9aU = 2c g g 3 J_ 
M A ' M A 



(A14) 

where we have defined 

u A = k(A\F^F^\A) -J2^ m ^ k ( A15 ) 
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APPENDIX B: EFFECTIVE POTENTIAL FOR HIGGS-ULTRALIGHT-SCALAR 

MIXING 

As discussed in section VI, ordinary matter can couple to the ultralight scalar 0, which 
mediates a long range WEP violating force, via its mixing with the Higgs. Here we show 
the computation of the effective potential which generates this mixing after electroweak 
symmetry breaking. This effective potential is generated at one loop via the sum of diagrams 
shown in Fig. 7 for the scalar singlet x an d real scalar triplet x models discussed in section 
VI. Working in unitary gauge where H = h/y/2 and in d-dimensions, one can write the sum 
of all diagrams in Fig. 7 as 



-iv e Us,h) = -i^sJ^Y. 



(a 2 h 



2\n 



-ing x S 



(2Tl) d ^ (k 2 
v ' n=0 v 

d d k 



M 2 )«+i 



1 



^ 2 (M 2 + a 2 h 2 ) 



E {27r) d (k 2 + M 2 + a 2 h 2 ) 
1 



- 7^ + ln47T + 1 - In 



M 2 



a 2 h/ 



/// 



(Bl) 



where the first line is obtained after performing a Wick rotation to Euclidean momentum 
space. The superscript in Vjg(S, h) denotes that it is only the part of the effective potential 
linear in S. We see from the above result that the coefficient of the S and Sh 2 operators are 
UV divergent. These divergences are understood from the need to renormalize the tadpole 
graph of S and the renormalizable coupling 5± of Eq.(38), corresponding to the first two 
diagrapms in Fig. 7. The remaining diagrams mix into non-renormalizable operators and 
are finite. The counterterms needed to cancel the UV divergences are 



i6V eB (S,h) = S 



" 16tt 2 



-7£ + ln47r +6i(/x) 



+ Sh 2 



16tt 2 



7s + ln47r + 



where and are scheme dependent finite quantities. 
The quadratic terms in the potential is given by 

U quad = \^\h 2 + v 2 s S 2 + ri s hS), (B3) 

as first shown in Eq.(40). As seen from Eqs.(76) and (80) the mixing angle for Higgs- 
ultralight-scalar mixing is given by 

sin0~^T~% (B4) 



45 



and we can write 



As = 2 



d 2 V r 



quad 



d 2 V eS (h,S) 



dSdh 



dSdh 



h=v,S=0 



g x a 2 



Ml 



47T 2 \ fJ, 

where we have defined the renormalized effective potential V c g as 



a| 9x v3 
4tt M2 : 



'eff- 



(B5) 



(B6) 



The first term with square brackets in Eq.(B5) corresponds to the renormalized value of 
v5i and the last term corresponds to the finite contribution from all non-renormalizable 
operators. This can be compared to Eq.(76) where we have included only the contribution 
from the renormalized 5\ coupling and the third diagram in Fig. 7 whose Wilson coefficent 
is denoted as C 2 . The above result, which is given in Eq.(79) of the text, is the generalized 
result where the contribution of the entire tower of higher dimension operators is resummed. 
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